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In the recent experiments like DIS at HERA or the heavy-ion experiments at RHIC, and also in 
expected LHC at CERN, the number of involved partons is very large, due to the high energy 
and/or the high number of participants of those experiments. These high parton densities 
should in principal lead to an extremely huge multiparticle production, but experimentally we 
have seen that this is not the case. So there should be a mechanism that reduces the number of 
created particle. Here, I review the problem of parton saturation and its implications through 
three in principal different approaches, but somewhat related: saturation in a geometrical 
approach, QCD saturation through the Color Glass Condensate and perturbative Pomeron 
approach with initial conditions. 

1 Geometrical approach to saturation: String models and percolation 

In many models of hadronic and nuclear collisions, color strings are exchanged between the 
projectile and the target. Those strings act as color sources of particles through the creation of 
q — q pairs from the sea. The number of strings grows with the energy and with the number of 
nucleons of the participant nuclei. 

In impact parameter space these strings are seen as circles inside the total collision area. 
When the density of strings becomes high the string color fields begin to overlap and eventually 
individual strings may fuse, forming a new string -cluster- which has a higher color charge 
at its ends, corresponding to the summation of the color charges located at the ends of the 
original strings. The new string clusters break into hadrons according to their higher color. 
As a result, there is a reduction of the total multiplicity. Also, as the energy- momenta of the 
original strings are summed to obtain the energy-momentum of the resulting cluster, the mean 
transverse momentum of the particles created by those clusters is increased compared to the one 
of the particles created from individual sources. 

As the number of strings increases, more strings overlap. Some years ago, it has been 
proposed in Ref. ^ that above a critical density of strings percolation occurs, so that paths of 



overlapping circles are formed through the whole collision area. Along these paths the medium 
behaves like a color conductor. Also in^^, we have made the remark that several fused strings 
can be considered as a domain of a non thermalized Quark Gluon Plasma. The percolation gives 
rise to the formation of a non thermalized Quark Gluon Plasma on a nuclear scale. 

Note that here we are not speaking about a final state interaction phenomenon, since there 
is no thermalization involved. In fact, what we are trying is to determine under what conditions 
the initial state configurations can lead to color connection, and more specifically, if variations 
of the initial state can lead to a transition from disconnected to connected color clusters. The 
results of such a study of the pre-equilibrium state in nuclear collisions do not depend on the 
subsequent evolution and thus in particular not require any kind of thermalization. 

The formalism underlying the transition from disconnected to connected systems is given by 
percolation theory, which describes geometric critical behavior. 

Consider placing small circular discs (color sources, strings or partons) of radius r onto a 
large circular manifold (the transverse nuclear plane) of radius R; the small discs may overlap. 
With increasing density 

NTrr^ 
^=^' 

this overlap will lead to more and larger connected clusters. The most interesting feature of 
this phenomenon is that the average cluster size increases very suddenly from very small to very 
large values. This suggests some kind of geometric critical behavior. In fact, the cluster size 
diverges at a critical threshold value rjc of the density. This appearance of an infinite cluster 
at rj = rjc is defined as percolation: the size of the cluster reaches the size of the system, rjc 
has been computed using Monte Carlo simulation, direct connectedness expansion and other 
different methods. All the results are in the range rjc = 1.12 1.175. 

In our model^ we had proposed a fixed radius for the independent color sources of r = 
0.2 0.25 fm, that corresponds to a momentum around 1 GeV. This value has been obtained 
from Monte Carlo simulations in the framework of the String Fusion Model Code (SFMC)I^ 
made at SPS energies. According to eq. (0), in order to estimate the density r/, one needs to 
know the number of sources A^. In our model, it is obtained from the SFMC, that, for nucleus- 
nucleus collisions, takes into account two contributions: one proportional to the number of 
participant nucleons -valence-like contribution- and another one proportional to the number of 
inelastic nucleon-nucleon collisions. Note that A^ will depend on the energy ^/s (or equivalently, 
on x) and on the number of participant nucleons A, so in some way the condition to achieve 
percolation depends on A and s, rj = r]{A,x). vri?^ corresponds simply to the nuclear overlap 
area, Sa_, at the given impact parameter. This overlap area can be determined in a Glauber 
study, using Woods-Saxon nuclear profiles. That leads to the following results: In our model, at 
SPS energies, the critical threshold for percolation could eventually been achieved for the most 
central Pb-Pb collisions, and for sure in Au-Au central collisions at RHIC energies and even in 
p-p collisions at LHC energies. Just as an example, in first approximation, one can estimates 
analytically that at very high energies, for central A- A collisions, if we take the number of 
initially created sources as proportional to the number of collisions, A^ oc A^^^ and the nuclear 
overlap area Sa = t^R^ oc A^^^ then 

^ ^ NsourcesSl ^ ^2/3 _ (2) 

2 The Color Glass Condensate 



Now we arrive to another approach, the QCD saturation through the formation of a Color Glass 
Condensate (CGC) The idea is the following: At high energy, the QCD cross-sections are 
controlled by small longitudinal momentum gluons in the hadron wave function, whose density 



grows rapidly with increasing energy or decreasing x, due to the enhancement of radiative pro- 
cess. If one apphes perturbation theory to this regime, one finds that, by resumming dominant 
radiative corrections at high energy, the BFKL equation leads to a gluon density that grows 
like a power of s and in consequence to a cross-section that violates the Froissart bound. Nev- 
ertheless, the use of perturbation theory to high-energy problems is not obvious. In fact, the 
BFKL and DGLAP equations are linear equations that neglet the interaction among the gluons. 
With increasing energy, recombination effects -that are non-linear- favored by the high density 
of gluons should become more important and lead to an eventual saturation of parton densities. 

These effects become important when the interaction probability for the gluons becomes of 
order one. Taking ^^r^ as the transverse size of the gluon and ^^'^^ ^ as the density of gluons, 
the interaction probability is expressed by 

a-sNc xG{x,Q'^) . 

Equivalently, for a given energy, saturation occurs for those gluons having a sufficiently large 
transverse size ~ larger than some critical value l/Qs{x,A). So the phenomenon 

of saturation introduces a characteristic momentum scale, the saturation momentum Qs{x,A), 
which is a measure of the density of the saturated gluons, and grows rapidly with 1/x and A 
(the atomic number). The probability of interaction -that can be understood as "overlapping" 
of the gluons in the transverse space- becomes of order one for those gluons with momenta 
^ Qs{x-,A) where 

Ql{x, A) = asN. ^-^^^ = ^^"^g^^' . (4) 
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For ^ Q1{x-, A), the non-linear effects are essential, since they are expected to soften the 
growth of the gluon distribution with r = ln(l/j;). For a nucleus, xGa{x, Q1) oc A and vri?^ oc 
A^/^ , so eq. @ predicts oc A^/^ . One can estimate the saturation scale by inserting the BFKL 
approximation into eq. @. This gives (with 5 ~ 1/3 and A ~ ca^ in a first approximation): 

Ql{x,A) ~ A'x-\ (5) 

which indicates that an efficient way to create a high-density environment is to combine large 
nuclei with moderately small values of x, as it is done at RHIC. In fact the estimated momentum 
for saturation at RHIC will be = 1 2 GeV, in accordance with the result of the previous 
section. 



3 Perturbative QCD pomeron w^ith saturation in the initial conditions 

Let us now try a different approach^!. Consider now the nucleus-nucleus interaction as governed 
by the exchange of pomerons. Its propagation is governed by the BFKL equation. Its interaction 
is realized by the triple pomeron vertex. Equations which describe nucleus-nucleus interaction 
in the perturbative QCD framework have been obtained in Knowing that the AGK rules are 
satisfied for the diagrams with BFKL pomerons interacting via the triple pomeron vertex one 
can conclude that the inclusive cross-section will be given by the convolution of two sums of fan 
diagrams propagating from the emitted particle towards the two nuclei. 

Taking A = B and constant nuclear density for \h\ < R^, one can find the inclusive cross- 
section in perturbative QCD as 

I^{y,k)=A^/\Rl^^ I d\e"^^[A^A{Y -y,r)][A<^>Aiy,r)], (6) 



where A is the two-dimensional Laplacian and ^(y, r) is the sum of all fan diagrams connecting 
the pomeron at rapidity y and of the transverse dimension r with the colliding nuclei, one at 
rest and the other at rapidity Y. The function (/>yi(y,r) = r)/(27rr^), in the momentum 
space, satisfies the well-known non-linear Balitsky-Kovchegov equation 

^ = -ff«„,,)-*^fe,,). (7) 

where y = ay, a = agNc/ir, as and Nc are the strong coupling constant and the number of 
colors, respectively, and H is the BFKL Hamiltonian. This equation has to be solved with the 
initial condition at y = determined by the color dipole distribution in the nucleon smeared by 
the profile function of the nucleus. 

We can take the initial condition in accordance with the Golec-Biernat distribution that 
takes into account saturation according to the CGC. Then we can find the following result: One 
can observe that whereas at relatively small momenta the inclusive cross-sections are propor- 
tional to A, that is to the number of participants, at larger momenta they grow with A faster, 
however noticeably slowlier than the number of collisions, approximately as A^'^. The interval 
of momenta for which 1^ (x A is growing with energy, so that one may conjecture that at infinite 
energies all the spectrum will be proportional to A. 

Note that in the string models, like Quark Gluon String Model or Dual Parton Model^each 
string corresponds to the exchange of pomerons, and the interaction among the strings would 
correspond to interaction among the pomerons through the triple pomeron vertex. 



4 Conclusions 



We have compared different models that takes into account saturation in different ways: from the 
semi-phenomenological fusing color sources picture for the soft domain including percolation, 
the QCD saturation through the Color Glass Condensate and those which follow from the 
pomeron approach, perturbatively derived from QCD, taking into account saturation in the 
initial conditions. In fact, it seems that the exchanged of elemental objects, color sources - 
strings, partons or pomerons-, should lead to a saturation in the initial conditions when the 
densities are high enough. 
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